Scattering theory of cooling and heating in opto-mechanical systems 
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We present a one-dimensional scattering theory which enables us to describe a wealth of effects 
arising from the coupling of the motional degree of freedom of scatterers to the electromagnetic 
field. Multiple scattering to all orders is taken into account. The theory is applied to describe the 
scheme of a Fabry-Perot resonator with one of its mirrors moving. The friction force, as well as the 
diffusion, acting on the moving mirror is derived. In the limit of a small reflection coefficient, the 
same model provides for the description of the mechanical effect of light on an atom moving in front 
of a mirror. 

PACS numbers: 37.10.De, 37.10.Vz, 42.70.Qs 



I. INTRODUCTION 



The use of light forces to manipulate mechanical mo- 
tion has been extended by now from the translational 
motion of single atoms [l|, 0] to the motional modes 
of massive systems, such as the oscillations of a micro- 
mechanical mirror [3, 0, 0]- The theoretical approach 
to describe the mechanical effect of light on the center- 
of-mass motion of atoms is completely distinct from the 
models dealing with vibrating optical resonators. In the 
first case, theories are based on the assumption that 
atoms are very weak scatterers in free space, negligibly 
perturbing the impinging bright laser beams @. In the 
other case, the influence of the moving massive compo- 
nent on the radiation field is so strong that it is con- 
sidered a (moving) boundary condition defining a single 
or a few modes of the field participating in the opto- 
mechanical coupling 0, [H . This is clearly the case for a 
Fabry-Perot type resonator with one of its mirrors mov- 
ing [g, [l(| OjJ. We argue that these two cases can be 
dealt with as two extremes of a general system that can 
be described in a unified theoretical framework. 

In this paper we develop and present a scattering 
theory for opto-mechanically coupled systems, allowing 
for the efficient description of the motion of arbitrary 
combinations of atoms and mirrors interacting through 
the radiation field. We will restrict the model to one- 
dimensional motion and small velocities. The main build- 
ing block is the beamsplitter transfer matrix (l^.[l3|. i.e., 
the local relation between light field amplitudes at the 
two sides of a scatterer. We will calculate the radiation 
force acting on a moving scatterer up to linear order in 
the velocity. The model is completed by including the 
quantum fluctuations of the radiation force which stem 



from the quantized nature of the field. We will determine 
the momentum diffusion coefficient corresponding to the 
minimum quantum noise level. 

The system we will consider in some detail is com- 
posed of two mirrors; one of them is fixed in space, whilst 
the other one is mobile. This is the generic scheme for 
radiation-pressure cooling of moving mirrors [hi . [Tol fioj ] . 
At the same time, in the limit of low reflection the moving 
mirror can equally well represent a single atomic dipole 
interacting with its mirror image in front of a highly re- 
flecting surface [13, d, Ell- The scattering model de- 
scription of this example gives a clear recipe for general- 
izing the method to more complex systems. 



II. MODEL 
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FIG. 1: The four different modes that interact through a 
point-like beamsplitter in ID. 

Consider a point-like scatterer (or beamsplitter), BS, 
moving along the s x' axis on the trajectory XBs(t)- Out- 
side the scatterer, the electric field E can be expressed 
in terms of a discrete sum of left- and right-propagating 
plane wave modes with different wave numbers, k, and 
hence different frequencies, to = kc: 



E* [Mk)e 
£ fe [C{k)t 



ihx—iwt 
— ikx—iujt 



B{k)e lk 
■D(k} 



ikx — iujtl 



Akx—iujtl 



CC. 
- C.C. 



*To whom all correspondence should be addressed. Electronic ad- 
dress: andre.xuereb@soton.ac.uk 
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where A(k) and B{k) are the mode amplitudes on the left 
side, x < XBs(t), while C(k) and D(k) are the amplitudes 
on the right side, x > xbs(£)j of BS. This is a simplifying 
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assumption and all our results also hold for a continuum 
of field modes. In accordance, the magnetic field is [2fJ] 
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As depicted schematically in Fig. [IJ the scatterer mixes 
these waves. Our first goal is the derivation of the trans- 
verse matrix M connecting the field amplitudes on the 
right to those on the left side of a beamsplitter moving 
at a fixed velocity v. This relation is well-known [l2j for 
an immobile scatterer. Therefore, let us first transform 
the electromagnetic field into a frame moving with the 
instantaneous velocity v of the BS. 



A. Transfer matrix for an immobile beamsplitter 

In the frame co-moving with BS, the interaction of the 
field with the scatterer at x' = can be characterized by 
the single parameter £ by means of the one dimensional 
wave equation [lj], H(| , 



The electric field can be considered in a modal decompo- 
sition similar to Eq. |T]). Since a fixed beamsplitter cou- 
ples only the plane waves with identical frequency and 
wave number, the stationary scattering can be fully de- 
scribed within the closed set of modes 



E'{x',t') = 
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B'e ikx '- iut ' + c.c. x' <0 
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+ c.c. x 1 > 0. 



where the index k has been dropped. A linear relation 
between the field amplitudes on the right of the scat- 
terer and those on the left can be derived from the wave 
equation, 
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In the second form of the transfer matrix Mo, we ex- 
pressed it in terms of the reflectivity r and transmissivity 
t of the beamsplitter. This latter form is more conve- 
nient to describe moving mirrors, while for atoms the 
scattering strength parameter C can be readily expressed 
in terms of the polarizability [12j, 
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e Q XS ' 



where a is the linear polarizability and S is the effective 
beam cross section. For a two-level, unsaturated atom 



with transition frequency u>a and linewidth T (HWHM), 
for example, 
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where cta = 4p- is the resonant radiative cross section 
of an atom. In this case the transfer matrix depends on 
the wave number k, which might lead to essential effects, 
e.g., Doppler cooling, close to resonance with the atom 
(see Section QTE}. 



B. Transfer matrix for a moving beamsplitter 

The transformation back into the laboratory-fixed 
frame involves the change of the coordinates, x' = x — vt 
and t' = t, and the Lorentz-boost of the electric field up 
to linear order in v/c [20, §11.10]: 

E = E' + vB' , 

where we assumed that E and E are polarized in the 
'j/' direction, B and B are polarized in the l z' direction, 
and the velocity is along the x axis. The electric field in 
the laboratory frame becomes 



(x— vt) — iuj't 



E(x,t) = Y,\ A'{k')e- ik ' (x - vt) - iuJ ' t +B'{k')e lk ' 
y { 

-- A'{k')e- lk ' {x - vt) - lu ' t ^B'{k')e lk ' { - x - vt) - w ' t j+c.c. 
= ^ (l — ~) A' (k + kv/c) e -^a+v/c)x-zut 

k 

+ (1 + f ) B' (k - kv/c) e iW-"/c)x-iwt ^ 

which can be expressed as a linear transformation L(v) 
of the amplitudes, 



\B{k))- L ^ \B'{k)) 



with 
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This construction is explored further in Appendix [Al 
Here we defined the operator P v : f(k) i— > / (fc + fe-), 
which represents the Doppler-shift of the plane waves in 
a moving frame. Obviously, = L(—v) to first 

order in -. The total action of the moving BS, 



C(k)\ ~ (A(k)\ 
D(k)J - {B(k)J 



can then be obtained from 
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where we have assumed that r and t do not depend on 
the wave number. 

Compared to Mo in Eq. ((4}, the difference lies in the 
off-diagonal terms including the Dopplcr-shift imposed 
by the reflection on a moving mirror. In other words, 
the coupled counter-propagating plane wave modes differ 
in wave numbers, i.e., k (l + -) right-propagating waves 
couple to —k (l — |) left-propagating waves. Further- 
more, if the polarizability itself depends on the wave 
number k, e.g., as in Eq. (0, the Doppler-shift opera- 
tor acts also on it. To make this effect explicit, to linear 
order in -, M can be written as 
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The transfer matrix in the laboratory frame can thus 
be conceived as a 2-by-2 supermatrix acting also in the 
/c-space. The amplitude C at a given wave number k, 
i.e., C(k), is combined with the amplitudes A(k) and 
B{k — 2fc~). A similar statement holds for D(k). 

Starting from the knowledge of the incoming field am- 
plitudes, this transfer matrix allows for calculating the 
total electromagnetic field around a beamsplitter mov- 
ing with a fixed velocity. In the next step, we derive the 
force on the moving scatterer through the Maxwell stress 
tensor. 



C. Force on a medium in an electromagnetic field 

The Maxwell stress tensor (see 20, §6.7]) is defined, 
for a homogeneous medium in one dimension, x, as 



T = -^\E\ 



c \B\ 



where the electric field E and the magnetic field B, 
Eq. |T]) and Eq. ([2|), respectively, have no components 
along x. It is trivial, then, to see that after applying the 
rotating wave approximation, we obtain 

T xx = -2e | A(k)e- lkx ~ iut + 1 ^ B{k)e ikx ~ lL 

. k k 

since the cross terms in \E\ 2 and \B\ 2 have opposite signs. 
Note that T varies on time scales of the order of the op- 
tical period. Let us now introduce a characteristic time, 
t, over which the variations in T will be averaged. At 

x = 0, 
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In the approximation we assumed that the frequency 
bandwidth of the excited modes, A = max {tUi — loj}, 
around the central frequency, ujq , is so narrow that At <C 
2ir. Since the broadening is due to the Doppler-shift, 
A ~ 2a>o-j, where v is the speed of the beamsplitter. For 
example, taking v to be the typical speed of atoms in a 
magneto-optical trap, we require r <C 7r/ (ljo-j ~ s. 
The time needed to reach the stationary regime of scat- 
tering is typically much shorter and thus this condition 
imposed on the averaging time r can be safely fulfilled. 

The force on the medium is given by the surface inte- 
gral of T on the surface, S, of a fictitious volume V = S 51 
enclosing the medium, where S is the mode area and 61 
the infinitesimal length of the volume along the 'x' axis. 
Then, this force is given by 



F = d) T xx n x dS 

= S[T xx (x^0+)-T xx {x^0-)], (8) 

where n x = sgn(ir) is the normal to S. Substituting 
the relevant expressions for T into the preceding formula 
gives 
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where A = [hw / (2Se c)]^ 1/2 J2 k A ( k ) is the photo- 
current amplitude, and similarly for B, C and D, their 
modulus square giving the number of photons crossing a 
unit surface per unit time. Although we considered first 
the electric field composed of independent modes, in the 
force expression only the sums of the mode amplitudes 



D. Quantum fluctuations of the force 

In the previous subsection the force was derived based 
on the assumption that the field amplitudes are c- 
numbers. In order to describe the inherent quantum fluc- 
tuations of the force, we need to resort to the quantum 
theory of fields and represent the mode amplitudes by 
operators: A(k) — > A(k). To leading order the fluctu- 
ations of the force acting on a beamsplitter amount to 
a momentum diffusion process [2ll I22I ]. The diffusion 
coefficient will be evaluated in the following in the min- 
imum, quantum noise limit, which occurs in the case of 
coherent-state fields [231 ] . 

The diffusion coefficient can be deduced from the 
second-order correlation function of the force opera- 
tor MM 



(F(t)F(t>))-(F(t)) 2 = D(t)6(t-t'). 



(10) 



The evaluation of this quantum correlation is system spe- 
cific. Quantum correlations, i.e., the operator algebra of 
the mode amplitudes A{k), B(k), C(k), and D(k), are 
influenced by multiple scattering and thus depend on the 
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total transfer matrix of the entire system. The simplest 
case is a single beam splitter where the "input" modes 
B(k) and C(k) have independent fluctuations. The cal- 
culation, delegated to Appendix [Bj includes all the steps 
needed for the treatment of a general system. The diffu- 
sion coefficient for a single beamsplitter is obtained as 



D = (hk) 2 (\ 



A\ +\B\ +\C\ +\D\ 
+ 2Re{rA*B ~iA*C} 

+ 2Rc{xD*C -tD*B}\ , (11) 



where A, B, C, D are the photo-current amplitudes (their 
modulus square is of the units of 1/sec), obeying Eq. ([6j) 
for v — 0. 

As an example, let us consider the diffusion coefficient 
for a two-level atom illuminated by counter-propagating 
monochromatic light waves. Using the polarizability £, 
the transmission and reflection coefficients can be ex- 
pressed as t = 1/(1 — i() and r = — «C)> respectively 
(see Eq. (j4j). Eq. (fTTj) can then be rewritten in the form 
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(12) 

where the first term, apart from the factor |1 — i(\ 2 , cor- 
responds to the result well-known from laser cooling the- 
ory, as shown in the next section. Note that the diffu- 
sion process due to the recoil accompanying the sponta- 
neous emission of a photon (see |25j |) is missing from this 
result — the detailed modeling of absorption, i.e., scat- 
tering photons into the three dimensional space, is not 
included in our approach. 



E. Example: Force on a moving beamsplitter 

We will now use Eq. © to derive a general expres- 
sion for the force on a moving beamsplitter illuminated 
by two counterpropagating, monochromatic, plane waves 
with amplitudes Bq and Co- On using Eq. © to ex- 
press the outgoing field modes in terms of the incoming 
ones, we note that the outgoing amplitudes comprise two 
monochromatic terms each: 



1(2 
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B + Co 
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These relations are substituted into Eq. ([9]), giving 

^ = {2^/|i-c(i + f|§)| 2 } 

x|(lm{C}+|C| 2 + i>^ ! )(l^o| 2 -|Co| 
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+ 2 (| W Im{c*§} - Re{C}) Im{B C *} 
(2Im{C} 



Re{B C*} 



(13) 



For v = this result reduces to the one in [1 31 ] - 
Most of the w-dependent terms arise from the frequency- 
dependence of the polarisability. These are the dominant 
terms in the case of a quasi- resonant excitation of a res- 
onant scatterer, such as a two-level atom, since the pref- 



actor 



~C~5u> '~ T ex P resses resonant enhancement. The 
w-dependent terms linear in the polarizability £ are in 
perfect agreement with the friction forces known from 
standard laser cooling theory, both for propagating and 
for standing waves. For example, assuming identical laser 
powers from the two sides, giving a standing wave with 
wavenumber fco , and averaging spatially gives 



-4Kk%\B Im 



{di} 



(14) 



for small £ and to first order in ^, which can be immedi- 
ately recognized as the friction force in ordinary Doppler 
cooling [6( when one uses the definition of C in Eq. J5J. 
Finally, by making similar substitutions into Eq. (|12[) . we 
obtain 



D = 8(/ifco) 2 Im{C} \B \ sm 2 (k x) 



(15) 



and 



which, excluding the diffusion effects due to spontaneous 
emission, matches the standard result in [25j. Note, 
however, that the scattering theory leads to a more 
general result which is represented by the terms of higher 
order in (. These terms describe the back-action of the 
scatterer on the field, an effect neglected in free-space 
laser cooling theory. 

The general result in Eq. (flU)) reveals that this velocity- 
dependent force also acts on a scatterer whose polaris- 
ability is independent of the frequency. This is a very 
general class and we will only focus on such scatterers in 
the following. 



III. GENERAL SYSTEM OF A FIXED AND A 
MOBILE SCATTERER 
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Consider the model in Fig. [2] where the scatterer, or 
'atom', has a polarizability Q uniform over the frequency 
range of interest. Letting M a , M p and M m be the transfer 
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FIG. 2: Physical parameters of our model. A, B, etc. repre- 
sent the field mode amplitudes. 



matrices for the atom, propagation and mirror, respec- 
tively, we obtain the relation: 



A4(fc)\ fC(k) 
I ) ~ MaM P M - I D(k) 



M a = 



, where 
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, and M m = - 



t 4 - r 
— r 



The distance between the atom and the mirror is denoted 
by d. Note that the free propagation transfer matrix M p 
is non-uniform in the fc-space, and therefore the Doppler- 
shift has an influence on the phase shift accumulated be- 
tween two scattering events. 

The boundary condition is set as follows. Since there 
is no incoming field from the right, C(k) = for all 
k. The incoming field from the left is assumed to be 
monochromatic, B(k) = B5(k — ko), with kg being the 
pump wavcnumber. The resulting field comprises modes 
with wavenumbers in a narrow region around ko. In 
the laboratory frame the field mode A(k) interacts with 
B(k - 2fc|) and C'(k) through the Doppler-shift, and 
similarly for D'(k). From C(k) = it directly follows 
that 
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(16) 



We will need the sum of amplitudes, A = j A(k) dk/B, 
defined relative to the incoming amplitude B = 
/ B{k)dk. Note that / P v f(k) dk = J f(k) dk. Thus, 
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It is worth introducing the reference point at a distance 
L = 2Nir/ko from the fixed mirror, where the integer 
N is such that the moving atom's position x is within 
a wavelength of this reference point. Then the atom- 
mirror distance can be replaced by d = L — x, and k$L 
drops from all the trigonometric functions. The solution, 
Eq. (|17p . has a clear physical meaning, in that the re- 
flected field, A, can be decomposed into an interfering 
sum of fields: the first term is the reflection directly from 
the atom, whereas the summation is over the electric field 
undergoing successive atom-mirror round-trips. We can 
also write the preceding expression in closed form: 
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-2ikox 
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-2ikQX 
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(l - i( - xiCe~ 2ik ° x ) A 



(18) 



This result is valid for arbitrary Q. The main virtue of our 
approach is clearly seen, in that we can smoothly move 
from C = 0, which indicates the absence of the mobile 
scatterer, to |C| — > oo, which corresponds to a perfectly 
reflecting mirror, i.e., a moving boundary condition. 

Let us outline some of the generic features of the above 
calculation that would be encountered in a general config- 
uration of scatterers. By using the formal Doppler-shift 
operators, we benefit from the transfer matrix method in 
keeping the description of the system as a whole within 
two-by-two matrices. The input-output relation for the 
total system is always obtained in a form similar to that 
of Eq. (flTjj) . As long as the Doppler broadening is well 
below the transient time broadening of the system, the 
calculation of forces and diffusion requires solely the sum 
of the mode amplitudes. An important point is that 
the integrated action of the Doppler-shift operator P v 
on monochromatic fields is a shift in fc-space. Therefore, 
by interchanging the order of terms and putting the P v 
terms just to the left of the input field amplitudes, they 
can be eliminated, such as in Eq. (|17[) . Finally, up to 
first order in v/c, the resulting power series, a trace of 
multiple reflections, can be evaluated in a closed form, as 
shown in Eq. (p~8|) . In conclusion, the illustrated method 
lends itself for the description of more complex schemes, 
for example, the cooling of a moving, partially reflective 
mirror in a high-finesse Fabry-Perot resonator [26| . 
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A. Force acting on the mobile scatterer 

To obtain the force on the moving scatterer, we also 
need to evaluate C'(k) and D'(k): 



C'{k) 
D'(k) 



1-iC 

i( (1 + 2^) P.' 1 



iC(l-2f)fW 

i + »C 



(A(k)\ 
{B(k)J 

(19) 



where we applied the inverse of the transfer matrix M a . 
Next, we make the following definitions: 



c = - 



B = \B\ 2 , 

2 

C'(k) dk , 



D'(k) dk 



and a simple calculation leads to 



CHl-zC| 2 A+K(l-2^)| 2 

+ 2He{»C*(l-*C)(l-2^M}, 

D = |<(l + 2f)| 2 A+|l + iC| 2 

+ 2Rc{i((l + i(*)(l + 2^)A} . 

Thereby the force acting on the scatterer is obtained as 

F = (hu/c)M (A + 1 - C - D) 

= -2ftfc l( [|C| 2 (1 + 2f ) + Im{C}] A 
+ |C| 2 (l-2f)-Im{C} 
+ 2Re{i((l-i()A}) , (20) 

where A has to be substituted from Eq. (fP8|) . The coef- 
ficient of the term linear in velocity, the 'friction coeffi- 
cient' (3, is plotted in Fig. [3] as a function of the position x 
in a half- wavelength range for various values of £. When 
varying the coupling strength from £ = 0.01 up to Q = 1, 
the friction coefficient transforms between two charac- 
teristic regimes. For small coupling the linear velocity 
dependence tends to a simple sinusoidal function while, 
for large coupling, the friction exhibits a pronounced res- 
onance in a narrow range. This resonance arises from 
the increased number of reflections between the mobile 
scatterer and the fixed mirror. It can be observed that 
the resonance shifts towards k$x — n on increasing £. 
In the opposite limit of small £, the maximum friction 
is obtained periodically at (n — j)tt/2 according to the 
sinusoidal function. The position of the maximum fric- 
tion is plotted in Fig. 4(a) showing the transition from 
77r/8 to 7r. The maximum friction force is plotted in 
Fig. |4(b)[ showing the two limiting cases of £ 2 behavior, 
in the limit of small £, and £ 6 behavior, in the limit of 
large (. These two cases are described in section HVl and 
section [V] respectively. 
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FIG. 3: (Color online.) The position dependence of the lin- 
ear coefficient of the velocity-dependent force acting on the 
mobile scatterer in Fig. [2] for various scattering parameters 
£ evaluated by using Eq. (Tig)) and Eq. (J20jl with k L = 100. 
The fixed mirror is assumed to be a perfect mirror. In order 
to fit all the curves into the same range, they are divided by 
the factors indicated in the figure. 



B. Diffusion coefficient 

The calculation of the diffusion coefficient proceeds 
along the same lines as that corresponding to a single 
beamsplitter, shown in Appendix |E] The difference is 
that the modes B(k) and C'(k) around the mobile scat- 
terer are not independent, for the reflection at the fixed 
mirror mixes them. Therefore, all the modes A, B, C, 
and D' have to be expressed in terms of the leftmost and 
rightmost incoming modes, B(k) and C(k), respectively. 
Instead of the derivation of such a general result for the 
diffusion, here we will restrict ourselves to the special 
case of r = —1 (<=> perfect mirror) and real C (^ no ab- 
sorption in the moving mirror). In this special case the 
diffusion calculation simplifies a lot, because (i) the per- 
fect mirror prevents the modes C from penetrating into 
the interaction region, and (ii) quantum noise accompa- 
nying absorption does not intrude in the motion of the 
scatterer. 

Only the modes B(k) impart independent quantum 
fluctuations. When all the amplitudes around the scat- 
terer are expressed in terms of B(k), and are inserted 
into the force correlation function given in Eq. (fit)]) , the 
commutator [b(i), ^(t 1 )] appears in all the terms (see Ap- 
pendix |5| . Straightforward algebra leads to 



D = h 2 k£M(A + 1 - C • 



(21) 



We emphasize that the above result is not general: the 
diffusion is not necessarily proportional to the square of 
the force. This simple relation here follows from the as- 
sumptions, r = — 1 and Im{£} = 0, declared above. 

To be consistent with the calculation of the friction 
force linear in velocity, the diffusion should be evaluated 
only for v = 0. From the ratio of these two coefficients, 
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FIG. 4: (a) The position of the maximum friction force, kox m , 
as a function of the dimensionless scattering parameter £ (on 
a semilog scale) acting on the scatterer in Fig. [2] with the 
fixed mirror being a perfect mirror. This position shifts from 
7n/8 to 7r on increasing (b) A similar plot, showing the 
maximum friction force as a function of ( (on a log-log scale) 
with koL — 100. In the limit of small the force scales as £ 2 
(cf. Eq. (|25|l ; dashed line) whereas in the limit of large zeta it 
scales as ( 6 (cf. Eq. (gTJ); dotted line). 



the steady-state temperature can be deduced. The veloc- 
ity independent components of the modes obey the fol- 
lowing relations: A = 1 and C = W (all incoming power 
is reflected). Therefore the diffusion coefficent further 
simplifies, 



D = Ah z k, 



2 ;„2ti 



1 -iC + Ke' 



-2ikox I 



(22) 



In Fig. [5l the temperature ksT = D/(2/3), where (3 
is the friction coefficient, is plotted as a function of the 
scattering parameter £. The friction and the diffusion 
coefficients are taken at the position where the friction 
is maximum, as shown in Fig. 4(a) The two limits of 



small and large scattering parameter £ will be analysed 
in section IIVI and section [Vj respectively. 
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FIG. 5: Characteristic temperature for the two-scatterer sys- 
tem of Fig. [2] given by the ratio of the diffusion and friction 
coefficients in the points where the friction is maximum, as a 
function of the dimensionless scattering parameter £ on a log- 
log scale. Constant and 1/C 2 dependence can be read off in 
the limits of small and large £, respectively. The fixed mirror 
is a perfect mirror. 



IV. ATOM IN FRONT OF A PERFECT 
MIRROR 

An atom pumped with a far off-resonance beam can be 
modelled as a moving mirror with small and real In 
this section we accordingly truncate our expressions to 
second order in We also assume that the fixed mirror 
is perfect; i.e., t = —1 and t = 0. Thus, 



F = 2/ifc B{2CIm{.4} - 2C, 2 Re{A} 

)A-C 2 (i-f)} 



C 2 (i 



(23) 



To obtain F to second order in £, we need A to first 
order. Using Eq. (17]) and Eq. (23}, we obtain: 



A = — e~ 2ik ° x 
+ Cf[ 

and 



C(i - 2ie 
2i + 2ie~ uk ° x 



2iknx 



4fc (£ 



-4:ikox 



■ x)e 



-4ikox~\ 



(24) 



F = 4ftfcoB(Csin(2fc x) 

- C 2 { 2 sin 2 (k x) [4 cos 2 (k x) - l] } 
-C 2 f [4sin 2 (2A: x) 

-4fc (i-a;)siri(4A;oa;)]), (25) 

in agreement with In the far field (x ^> A), the dom- 
inant friction term in the preceding expression is the last 
term, which renders the sin(4fcoa;) position dependence 
shown in Fig. [3] for £ = 0.01. 

We are now in a position to derive the diffusion 
coefficient for this system. By substituting Eq. 
into Eq. (|2"Tj) and setting v — 0, we obtain 

d = 8(hk ) 2 ( 2 m . 
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This allows us to estimate the equilibrium temperature 
for such a system at a position of maximum friction: 



2r 



, where r = 2(L — x)/c, 



(26) 



which we note is identical in form to the Doppler tem- 
perature for a two-level atom undergoing free-space laser 
cooling [6], but where we have replaced the upper state 
lifetime, l/(2T), by the round-trip time delay between 
the atom and the mirror. Note that this temperature 
corresponds to the constant value presented in Fig. [5] for 
C < 0.1. 



V. OPTICAL RESONATOR WITH MOBILE 
MIRROR 

After the small polarizability case of the previous sec- 
tion, we will now consider the |£| — > oo limit. We again 
assume that the fixed mirror of the resonator is perfect, 
with r = — 1, and that C = 0. For simplicity, we assume 
that the moving mirror has a real polarizability. We ex- 
pand the field mode amplitudes as power series in v/c, 
such that A = Ao + ~A\ + . . . , and similarly for C. 

Let us first calculate the field in the resonator for v = 0. 
It follows from Eq. (19]) that 

p —2itp 



C' = (l- iQAo -i( = 



1 - i( + i(e~ 2i f ' 
with ip — kgd, which has a maximum at ipo obeying 



tan(2iy9o) 



In the limit of £ — ► oo, the resonance is Lorentzian: 

p -2itp 

C' - 

— r 

2i(l - iQ - <po) - ij^2 

with a width of 1/(4C 2 ). 

The perfect mirror reflects the total power incoming 
from the left, B. Moreover, for real £, there is no absorp- 
tion in the moving mirror, so the outgoing intensity has 
to be equal to the incoming one: A = 1. This is true if 
v = 0; for v =/= 0, the field can do work on the mirror. 
The expansion of the back reflected intensity to linear or- 
der in velocity reads A = 1 + 2^ Re{^4g^4i}. Extracting 
the velocity-dependent terms for the general form of the 
force in Eq. (|2"D|) , it reduces to 

Ft = f4frfc BCIm{A/ (1 + i£ - i(e 2itp )} , 
which, after some algebra, leads to 



Fi = 



■\lhklL- 



l3* 



(27) 



On substituting k = c/(4L£ 2 ), Ac = — c((p — <po)/L, 
r] 2 /(2k) = B, and G = c 2 k 2 /L 2 , the friction force renders 
that derived from the usual radiation pressure Hamilto- 
nian in Appendix [Cj 

Expressing the field modes interacting with the mobile 
mirror in terms of the input field mode and performing 
a calculation similar to that leading to Eq. (fTTj) readily 
gives 



D ^ A{hk ) 2 \C' Q \ 



(hk ) 2 M 



4C 4 



4F2) +{ip-ip y 



The resulting temperature thereby attains a minimum at 
4£ 2 (</? — <po) = 1, i.e., Ac = —K, in analogy with free- 
space Doppler cooling, at which point we have 



fc R T: 



-5£_-IfiK 

8C 2 L ~ 2hK - 



(28) 



Again, this asymptotic behavior is reflected in Fig. [5] for 
large £• We note the similarity of the preceding expres- 
sion with the temperature of an atom cooled in a cavity, 
in the good cavity limit We conjecture that this 

is due to the fact that both systems can be considered 
to involve the coupling of a laser with a system having 
a decay rate K. This result also holds for the case of an 
atom undergoing mirror-mediated cooling, as can be seen 
in Eq. ([26]) . 

It is also important to note that the above discussion 
only treats the effects of the light fields on the scatterer. 
As such, the temperature limit, Eq. (|28p . is intrinsic to 
the light forces, and the mechanical damping and heating 
processes present in a real, macroscopic mirror-cooling 
setup are not taken into account. In practice, these heat- 
ing processes may dominate over the heating induced by 
the quantum noise in the light field 28, 29]. In such cases, 
radiation pressure cooling is a possible means to lower the 
equilibrium temperature owing to the additional, optical 
damping process. 



VI. CONCLUSIONS 

We have presented a powerful extension of an existing 
theoretical framework to analyse the interaction between 
light and matter. The theory we presented is based on 
the transfer matrix method for dealing with the interac- 
tion between scatterers and a light field, and is therefore 
able to handle complex optical systems, made from sev- 
eral elements, with relative ease. Through the use of the 
Maxwell stress tensor one can calculate the force acting 
on any of the elements in the system. We have gener- 
alized the transfer matrix for slowly moving scatterers, 
thereby the corrections first-order in v/c can be calcu- 
lated for the electromagnetic field as well as for the ra- 
diation force acting on the scatterer. Furthermore, one 
can express this force in terms of the operators repre- 
senting the quantized field modes interacting with the 
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scatterer and consequently derive the momentum diffu- 
sion of the scatterer due to the quantum noise present in 
the fields. Our scattering theory can also transparently 
cover the whole range of interaction strengths, from the 
perturbative interaction between a weak standing wave 
and a single atom to the very strong (quasi boundary- 
condition) interaction between a pump light field and a 
Fabry-Perot cavity with a moving mirror. 

We also applied this framework to three different 
laser cooling configurations: optical molasses, mirror- 
mediated cooling and cooling of micromirrors. We de- 
rived the forces on an atom arising from its interaction 
with the light field, as well as an estimate for the equi- 
librium temperature an ensemble of atoms is expected to 
reach through this interaction. In the case of optical mo- 
lasses, which corresponds to the well-known Doppler tem- 
perature limit, the theory provides for additional force 
and diffusion terms related to the effect of the back-action 
of the atom on the radiation field. Although for single 
atoms in free space this back-action is feeble, it is re- 
sponsible for the modification of equilibrium properties 
and for collective effects in large optical lattices 
32| | . In the latter cases of a moving scatterer in front of 
a fixed mirror, our results are valid for arbitrary scatter- 
ing strength, i.e., spanning the parameter range from a 
single atom to high-reflectivity mirror. 
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APPENDIX A: THE DOPPLER-SHIFT 
OPERATOR 



number, written separately as: 

A B (k) = ai B(k + 2kf) , and 
A c (k) = a 2 C{k). 

We can therefore express A(k) as 

A(k) = ai B(k + 2fc§) + a 2 C(k) . 

Defining P v by P v : f(k) i— > f(k + k-), we have 

A(k) = P 2v a 1 B{k) + a 2 C(k) . 

A similar expression, involving P~ = P~ v , holds for 
D(k). These two operators can then be introduced 
into Eq. ((3|) as part of the Lorentz transformation, and 
thus into the transfer matrix for the moving scatterer, 
giving rise to the form shown in Eq. ([7]). The resulting 
transformation, for the transfer matrix M , of a scatterer 
moving with velocity v can be written as: 



(1 



(1 







)P- 



M 



(1 



)P- 











(1 



)P, 



to first order in ^, where the ordering of the elements of 
M is as described in the text. Note that this relation is 
general, in the sense that the elements of M can depend 
on k (see section Hi El) . 

For any finite v, P v is trivially a bounded operator, 
having unit norm. This property follows from the impor- 
tant relation J P™ f(k) dk = J f(k) dk, for any function 
f(k) and any integer m. 

This operation can be generalized to n — 2, 3 dimen- 
sions. We define a new operator by <Si(v) : /(k) h-> 
f(k + ki^ei), where e, is the unit vector along the ith 
coordinate axis, v is the velocity vector of the scatterer, 
and x = (xi, x 2l . . .) for any vector x. In particular, we 
have P v = Si(vei). Now, let L(v) be the 2n x 2n matrix 
operator: 



(1 



L )srV) 

(1 






f)^i(v) 







(1 + sl)4-i (v) 



A -< 

B 



S 



c 



FIG. 6: Reflection and transmission of a moving scatterer. 
B and C are the input field modes, and A is the output field 
mode. A further output field mode ('D') is not drawn because 
it is not relevant to our discussion in this section. 



Consider the situation in Fig. [6l in the laboratory 
frame, where S is a scatterer and suppose that B and 
C are known. A{k) has contributions arising from both 
B{k + 2k-) and C(k), where k is any arbitrary wave 



Then, the transfer matrix for the scatterer moving with 
velocity v is given by 

L(-v)ML(v), 

where M is the original transfer matrix for the scatterer, 
obtained in a manner such as that used to obtain Eq. ([4]) , 
for example. The ordering of the elements of M is such 
that it acts on the vector ( A\ (k) , B\ (k) , A 2 (k) , . . . J : 



/Ci(k)\ 
D 1 (k) 
C 2 (k) 



V 



= L(-v)ML(v) 



B 1 (k) 
A 2 (k) 

V J 
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with j4j(k) being the outgoing mode and -Bi(k) the in- 
coming mode along the ith axis in the negative half-space 
(assuming that the scatterer is at the origin); and Cj(k) 
the incoming mode and Z?i(k) the outgoing mode in the 
positive half-space. 



APPENDIX B: QUANTUM CORRELATION 
FUNCTION OF THE FORCE OPERATOR 



In quantum theory, we need to replace the mode am- 
plitudes A{k) by operators A(k), and similarly for the 
B, C, and D modes. The cross correlation of these op- 
erators is not trivial because of the boundary condition 
connecting the mode amplitudes A(k), B(k), C(k) and 
D(k). The input modes C{k) and B(k) can be consid- 
ered independent, and the commutator is non- vanishing 
for the creation and annihilation operators of the same 
mode, e.g., 



B(k),B^(k') = C(k),C\k') 
B(k),C^(k')] =0, 



Sk.k' 



assuming a discrete mode index of k 7 and a quantisation 
volume V = SI with S being the mode area and I a 
fictitious total length of the space in one dimension. 

We consider only the v = case, since our expressions 
are accurate up to first order in v/c. In the quantum 
description, the linear relation for the output modes is 

A(k) = tC(k) + tB(k) + ^E 
D{k) = rC(fc) + tB(k) + , 

where the transmission t = I/M22 = 1/(1 — iQ, and 
reflection r = M12/M22 = *C/(1 — *C)i as above. The 
fictitious amplitude E represents the quantum noise fed 
into the system by the absorption. For 7 = 1— (|r| 2 + |t| 2 ) , 
this noise ensures that the output modes obey the same 
commutation relations as the input ones, namely 



A(k),A\k') = D(k), D^(k') 
A(k),D\k')] =0. 



huj 

2^V 



However, the linear dependence implies that commuta- 
tors between input and output mode operators are 



A{k),B\k') 



= x 
= t 



B(k),B\k') 
C(k),C\k') 



and similar relations hold for the cross-commutators with 
D(k). 

The proper treatment of quantum fluctuations and the 
derivation of correlation functions require that the ex- 
plicit time dependence be considered. Let us introduce 



the time- varying operators 

a(i) = ^i(fc)e— 4 , 

k 

and similarly for b(t), c(t) and d(t). It follows that 



a(t),a*(t') 



-iu>{t-t') 



tluJ 

2cqcS 



S(t - 1') . 



Here we used that the non-excited, vacuum modes also 
contribute to force fluctuations. Therefore the Fourier- 
type summation extends to a broad frequency range and 
yields a S(t — t') on the much slower timescale of interest. 
A similar commutation relation applies to the operators 
b(t), c(t), and d(t). The cross-commutators can be de- 
rived directly from those concerning the modes, e.g., 



2eqcS 



S(t-t') 



The force operator is 

F = S[T ax (x -»■ 0+) - T xx (a 
as before, where 



0-)] 



T xx (x^0) = 



-2e 



-2e r 



(Bl) 

0" 
0+ 



is the quantized stress tensor. Assuming that the field 
is in a coherent state, in all normally ordered products, 
the mode amplitude operators can be replaced by the 
corresponding coherent state amplitudes, which are c- 
numbers: e.g., A(k) -> A(k) and A^(k) -> A*(k). The 
force operator in Eq. (|BL[) is normally ordered in this way; 
therefore coherent-state fields render, as a mean value of 
the quantum expressions, the force Eq. ([9]) derived from 
the classical theory based on the definition Eq. ([S]) . Non- 
trivial quantum effects arise from non-normally ordered 
products, such as the 4th-order product terms of the 
second order correlation function of the force Eq. (| L0[) . 
These terms can be evaluated straightforwardly by in- 
voking the above-derived commutators to rearrange the 
product into normal order. As an example, consider 

{a\t)a{t)a) (t')a{t')) = {at {t)a\t')a{t)a{t')) 

-{^(m(t'))^5(t- 



For radiation fields in coherent state, the first term is can- 
celed from the correlation function by the (a^> (t)d(t')} 2 
term. The coefficient of 5(t — t') in the second term is 
in normal order and can be replaced by c-numbers and 
then calculated identically as the force in section Hi El 

(aHt)a(t)) *\Y,A(k)\ 2 = J^\A\\ 

in terms of the photo-current intensity \A\ 2 . 

Assembling all similar contributions, originating from 
the non- vanishing commutators [b, tf], [c,c'], [d 7 d'], 
[a, fet], etc., one obtains Eq. (fTTj) presented in section lTl Dl 
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APPENDIX C: MIRROR COOLING VIA THE 
RADIATION PRESSURE COUPLING 
HAMILTONIAN 

We describe a generic opto-mechanical system com- 
posed of a single, damped-driven field mode coupled to 
the motion of a massive particle, whose Hamiltonian is 
given by [Hil 

H = huj c a)a + itvnitfe-™* - ae lut ) 

+ + V(x) + HGa'ax . 
2m 

where a and cfl are the annihilation and creation opera- 
tors of the mode, x and p are the position and momen- 
tum operators associated with the motion and we drop 
the carets to signify expectation values. The mode is 
driven by a coherent field with an effective amplitude r\ 
and frequency u>. This Hamiltonian describes, for exam- 
ple, the radiation pressure coupling of a moving mirror to 
the field in a Fabry-Perot resonator. In this case the cou- 
pling constant is G = lu c /L, rendering the cavity mode 
frequency detuning lo c x/ L provided the mirror is shifted 
by an amount x. Since the cavity mode is lossy with a 
photon escape rate of 2k, the total system is dissipative. 
Thereby, with a proper setting of the parameters, in par- 
ticular the cavity detuning Ac = lu — luc, the mirror mo- 
tion can be cooled. We will determine the corresponding 
friction force linear in velocity. 

In a frame rotating at frequency u, the Heisenberg 
equation of motion for the field mode amplitude reads 

a = [j(Ac — Gx) — k] a + r/ . 

where the noise term is omitted. We assume that the 
mirror moves along the trajectory x(t) w x + vt with 
fixed velocity v during the short time that is needed for 



the field mode to relax to its steady-state. The varia- 
tion of a arises from the explicit time dependence and 
from the motion of the mirror. A steady-state solution is 
sought in the form of a ~ a^(x) + va^'(x). On replac- 
ing this expansion into the above equation, and using the 
hydrodynamic derivative — > + vJ^, one obtains a 
hierarchy of equations of different orders of the velocity 
v. To zeroth order the adiabatic field is obtained as 



-i(Ac - Gx) + k ' 
The linear response of a to the mirror motion is then 

« (1) = ^ #« (0) 

z(A<7 — Gx) — k ox 

_ irjG 

[-i(A c -Gx) + n] 3 ' 

The force acting on the mirror derives from the defining 
equation p = jr[H,p] = —hGa)a. The force linear in 
velocity is 

Fi = -2vhGRe\a^a^\ = ^ , 

where we used x — without loss of generality. It can 
be seen that mirror cooling requires that Ac < 0, i.e., 
the cavity resonance frequency is above the pump fre- 
quency. In this case, for efficient excitation of the field in 
the resonator, the frequency of the pump photons is up- 
shifted at the expense of the mirror's kinetic energy. This 
cooling force has been derived in section [Vj as a limiting 
case of the more general scattering theory. To check the 
perfect agreement between the two results, the quantity 
corresponding to rj can be deduced from the total field 
energy in the resonator for an immobile mirror, which is 
Swca (0)t a (0) here. 
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